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ABSTRACT 


An  extension  of  a  previously  developed  approximate  theory,  the 
effective  dispersion  theory,  for  wave  propagation  in  a  periodically 
layered  elastic  body  is  given.  While  the  previous  theory  was  valid 
only  for  wave  propagation  normal  to  the  layering,  the  present  theory 
removes  this  restriction  in  the  case  of  antiplane  strain  (SH-wave) 
motion. 
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INTRODUCTION 

Delph  and  Herrmann  [1]  have  recently  proposed  a  conceptually  sim¬ 
ple,  yet  accurate,  approximate  theory  for  one-dimensional  wave  propa¬ 
gation  normal  to  the  layering  in  a  periodically  layered  elastic  com¬ 
posite.  This  theory,  called  the  effective  dispersion  theory,  has  been 
shown  to  model  quite  well  various  features  of  the  exact  dispersion 
spectrum  for  wave  propagation  normal  to  the  layering  in  a  periodically 
layered  unbounded  body.  In  particular,  the  effective  dispersion  theory 
contains  the  stopping  band  between  the  first  two  Brillouin  zones,  a 
feature  not  present  in  other  approximate  theories. 

This  theory  is  extended  in  the  present  work  to  two  dimensions  for 
antiplane  strain  motions  characterized  by  two  independent  wave  numbers. 
Thus  we  are  concerned  here  only  with  S  H  waves.  When  the  direction 
of  propagation  is  normal  to  the  layering,  the  extended  theory 
reduces  to  the  original,  one-dimensional  theory. 

The  extension  is  carried  out  in  a  simple  manner  by  including  ad¬ 
ditional  terms  in  the  strain  energy  density  function.  The  spectral 
matching  procedure  utilized  in  [1]  is  then  applied.  A  notable  feature 
of  the  matching  procedure  is  the  retention  of  the  stopping  band  between 
the  first  and  second  Brillouin  zones  which  was  present  in  the  one¬ 
dimensional  theory. 

On  the  basis  of  calculations  made  using  a  set  of  representative 
material  parameters,  it  appears  that  the  extended  two-dimensional  theory 
yields  excellent  results  for  values  of  wave  numbers  measured  parallel 
and  normal  to  the  layering  of  up  to  the  width  of  one  Brillouin  zone. 
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The  theory  also  appears  to  give  acceptable  results  over  a  smaller  range 
of  wave  numbers  in  the  second  Brillouin  zone. 

GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  layered  body,  B,  is  assumed  to  be  composed  of  periodically 
alternating  layers  of  homogeneous,  isotropic,  linearly  elastic  material 
(Fig.  1).  The  layers  have  thicknesses  2h  ;  2h *  ,  densities  p  ;  p' 
and  elastic  constants  x,u  ;  •  Any  two  adjacent  layers  comprise 

a  unit  cell,  and  this  unit  cell  is  repeated  periodically  along  the  y- 
axis.  The  body  is  taken  to  be  of  unbounded  extent  along  the  z-axis, 
but  bounded  in  planes  parallel  to  the  x-y  plane  by  the  curve  C(x,y)  , 
which  we  assume  to  be  sufficiently  smooth. 

For  antiplane  strain  motion,  the  displacement  component  parallel 
to  the  z-axis  is  taken  to  be  the  only  non-zero  displacement  component 
and  is  assumed  to  be  a  function  only  of  time  and  the  x  and  y  coordinates. 
We  now  assume  that  the  displacements  in  the  N-th  unit  cell  are  given  by 

w(x,y,t)  =  wQ(x,y,t)  |yN=0  +  yN<»(x,y,t)  |yN=0 

(1) 

w'(x,y,t)  =  wQ(x,y,t)|yN=0 

Here  w  and  w'  are  the  displacements  in  the  layers  with  unprimed  and 
primed  constants,  respectively,  y^  and  y^  are  local  coordinates 
having  their  origin  at  the  midplane  of  the  layers  with  unprimed  and 
primed  constants,  respectively.  As  noted  in  [1],  the  validity  of  equation 
(1)  depends  partly  on  the  assumption  that  the  layers  with  unprimed  con¬ 
stants  are  sufficiently  stiff  compared  to  the  layers  with  primed  constants 
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so  that  all  the  deformation  is  essentially  confined  to  the  softer 
layers. 

Following  [1],  we  now  postulate  a  specific  form  for  the  strain  and 
kinetic  energy  density  functions.  In  doing  so,  it  will  be  convenient, 
because  of  the  increased  complications  introduced  by  the  extra  spatial 
dimension, to  adopt  the  formalism  utilized  by  Achenbach,  Sun  and  Herrmann 
[2]  in  dealing  with  a  similar  problem.  We  thus  assume  that  the  strain 
energy  density  function  is  given  by 

W  =  Aijk£eijek£/2  +  Bijk£eijYk£/2  +  Ci  jk£Yi  jYlc£/2 

+  Di  jkfmn^ijk^n^2  +  ^i  jk£mn<ijk7£iw/2  (2) 

+  Fi jk£mnei jke-&mr/2  * 

Here  A.  ,F. ...  are  constants  to  be  determined  and  e..  ,  y.  .  , 

i  j  kc  i  j  k£mn  l  j  i  j 

<ijk  ,  and  e^mn  are  various  microstructural  strain  measures  intro¬ 
duced  by  Mindlin  [3].  In  our  case,  only  certain  of  these  are  taken  to 
be  non-zero.  These  are,  following  the  notation  of  [2] 

e32  =  e23  =  3ywo/2 
e13  =  e31  =  3xwo/2 
y23  *  3ywo  '  * 

(3) 

0223  =  3yywo  "  V 
<223  =  ay* 

IC1 23  =  K321  =  ax*/Z 
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where 


All  other  microstructural  strain  measures  not  defined  explicitly  in  (3) 
are  assumed  to  vanish.  Furthermore,  we  assume  for  the  sake  of  simplic¬ 
ity  that 


A2323  =  A2332  =  A3223  =  A3232 

A1 31 3  =  A3131  *  A1 331  =  A3131 

A1323  =  A2313  *  A1 332  =  A3213  =  A3123  =  A2331  =  A31 32  =  A3231  =0 

B2323  =  B3223 

(4) 

B1 323  =  B31 23  =  0 

°223123  =  D1 23223  =  °223321  =  °321 223  =  0 

D1 231 23  =  D321 321  =  D1 23321  =  °321 1 23 

E1 23223  =  E321 223  =  0 

The  kinetic  energy  density  is 

T  =  b-j  Wq/ 2  +  b2«ii2/2  (5) 


where  the  superposed  dot  indicates  partial  differentiation  with  respect 
to  time.  The  assumed  forms  for  the  strain  and  kinetic  energy  densities 
are  now  substituted  into  Hamilton's  principle,  which  we  write  as 
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(T- W)  dt  dV  + 


JN 


dt  =  0 


(6) 
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As  usual,  the  variations  6wQ  and  6iji  are  assumed  to  vanish  at 
t  =  t  ,  tj  .  Wg  is  the  work  done  by  external  tractions  on  the  bounding 
surface  of  the  body.  Again  following  [2],  we  assume  this  to  be  given  by 

<SWg  =  /s  CP(x,y)  <SwQ  +  M(x,y)  +  Q(x,y)  D6w0]  dS  (7) 

Here  D  is  the  normal  derivative  operator  on  the  bounding  curve  C(x,y) 

given  by  D  =  n  3/a  +  n  3/a  and  n  and  n  are,  respectively,  the  x 
x  x  y  y  x  y 

and  y  components  of  the  unit  normal  vector  to  C(x,y)  .  Body  forces  are 
assumed  to  vanish. 

The  derivation  of  the  governing  equations  and  the  associated  bound¬ 
ary  conditions  now  follows  that  presented  in  [2],  which  we  will  not 
repeat  in  detail.  Briefly,  we  make  use  of  equations  (2),  (3),  and  (4) 
to  define  the  following  generalized  stress  quantities 

T 23  H  aW/ae23  =  (A2323  +  B2323/2^  3yWo  “  B2323^2 
t13  5  3W/3e13  =  A13133xwo 

°23  H  3W/3y23  =  ^ B2323^2  +  C2323^  3ywo  "  C2323^ 

(8) 

P1 23  2  3W/3ki 23  =  D1 231 23 V 

U223  =  3W/3k223  =  (D223223  -  E223223/2)  3y*  +  E2232239yyWo/2 
x223  H  3W/36223  =  ^ E223223/2  '  F223223^  3y*  +  F2232233yywo 
Equations  (7)  and  (8),  after  several  applications  of  the  divergence 


theorem,  then  yield  the  equations  of  motion 


(9) 


-blwo  +  V23  +  3xT13  "  3yyx223  +  Sy°23  =  0 
“V  +  3yw223  +  Vl23  "  3yx223  +  *23  =  0 


with  the  associated  natural  boundary  conditions  on  the  surface  bounded 
by  C(x,y) 


<-Vl3  ‘  "y(T23  *  “23  '  sy*223>  '  C<1  ’  nx>  Vx  ‘  Vy  (sy"x 
+  (I  •  "y>  VyJ  nyx223  *  ny  C'VyVto  +  <>  *  ny>  V223] 

*  C'nxnyVy  *  (1  ‘  "p  Y'p  *223  +  F(x'x))  5wo  =  0 
^”nxpl 23  ”  ny  ^p223  ”  x223^  5tj3  -  0 

^~nyx223  +  D5uq  »  0 


3xV 


(10) 


Note  that  the  last  of  equations  (10)  implies  that  Q(x,y)  must  vanish 
on  surfaces  where  n^  =  0  .  This  is  somewhat  analogous  to  a  restriction 
present  in  the  effective  stiffness  theory  [2]. 


DETERMINATION  OF  COEFFICIENTS 

We  have  now  to  determine  the  various  constants  in  the  expressions 
for  the  strain  and  kinetic  energy  densities  given  by  equations  (2)  and 
(5).  This  is  accomplished,  as  in  [1],  by  matching  the  approximate 
dispersion  spectrum  for  antiplane  strain  wave  propagation  in  an  un¬ 
bounded  body  to  the  exact  result  given  in  [4].  To  this  end,  we  rewrite 
equations  (9)  as 


(11) 


•bl"o  +  al3yywo  +  a23yyyywo  ‘  a33y*  '  a45»y*  +  a73xxwo  =  0 


■b2*  +  a33ywo  +  Vyyywo  ’  a5*  ’  a6 V  +  Vxx*  =  0 


where,  with  the  aid  of  equation  (4),  we  have  defined 


al  =  A2323  +  B2323  +  C2323  a5  =  C2323 


a2  =  "F223223 

a3  =  B2323/2  +  C2323 

a4  =  E223223/2  '  F223223 


a6  =  E223223  '  F223223  '  °223223 

a  -A  < 

/  1 31 3 

a8  =  D1 231 23 


Equation  (11)  is  similar  to  the  governing  equations  derived  for  the  one¬ 
dimensional  case  in  [1],  but  contains  in  addition  terms  proportional  to 

Vo  and  v  • 

In  order  to  examine  the  case  of  harmonic  wave  propagation  in  an 
unbounded  body,  we  assume  solutions  of  the  form 


i(k  x  +  k  v  -  at)  i ( k  x  +  k  v  -  ut) 

wQ  =  Ae  T  ;  ip  =  Be  y 


Substitution  of  (13)  into  (11)  yields  a  set  of  two  linear,  homogeneous 
equations  in  A  and  B  ;  the  requirement  that  the  determinant  of  the 
matrix  of  coefficients  vanishes  gives  the  dispersion  equation  in  the 


auA  +  (6kA  +  Yk2  +  6  +  $k2)  w2  +  £k2  +  ak2  +  0kA  +  pk2  k2 
♦  ,kx  +  *  tkx  k5  ■  0 


The  parameters  a  ,  6  ,  y  ,  3  ,  £  and  9  are  various  combinations  of 
the  constants  a-p...,  b2  as  defined  in  [1].  Additionally  we  have 

$  =  ”(^7^2  3gb-j  )/T  ;  0  =  a^ay/T  ;  p  -  (a^ag  -  3g3y)/T  ; 

2  05) 

t  =  a7ag/T  ;  e  =  -a2ag/T  ;  T  =  a2ag  -  a4 


In  [1],  expressions  for  the  parameters  a-j  ,a2>a3,a4,a5,b-j ,  and  b2 
were  determined  by  matching  the  approximate  dispersion  spectrum  for  prop¬ 
agation  normal  to  the  layering  to  the  exact  solution,  [4],  by  matching 
the  approximate  mode  slope  at  the  end  of  the  first  Brillouin  zone  to  the 
exact  result,  and  by  imposing  the  requirement  that  the  frequencies  re¬ 
sulting  from  the  approximate  dispersion  relation  be  real-valued  for  all 
values  of  the  wave  number.  Since  the  governing  equations  are  homogeneous, 
the  value  of  ag  was  arbitrarily  set  to  unity.  We  retain  these  pre¬ 
viously-developed  values  in  the  present  work. 

In  order  to  determine  the  values  of  the  constants  a^  and  ag  , 
we  again  adopt  a  matching  procedure.  We  choose  to  match  first  the 

quantity  to  the  exact  value  obtained  from  [4],  which  we 

Kx  (o=kx=ky=0 

denote  as  s^  .  From  (14),  we  find  that 


“=kx=ky=0 


/-a/6  =  Sg 


(16) 


The  other  relation  is  obtained  by  requiring  the  imaginary  branch 
between  the  first  two  Brillouin  zones,  which  was  present  in  the  one¬ 
dimensional  theory,  to  be  preserved  for  values  of  kx  other  than  zero. 
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This  is  accomplished,  as  in  the  one-dimensional  theory,  by  requiring  the 
quantity  3u/3k^  to  vanish  along  the  end  of  the  first  Brillouin  zone, 
where  k^  has  the  value  k^  =  kg  .  From  (14),  we  have  that 


VS 


-ks[(28k|  +  y) 


m2  + 


(2ekg  +  p) 


k2  +  3kJ  +  20k2  +  ?] 


'B 


j[2ao)  +  Sk, 


+  ykg  +  6  +  $k2] 


07) 


Now  from  the  definitions  of  0  ,  y  ,  9  and  £  given  in  [1],  we  find 


28k£  +  y  =  0 


3kg  +  29  kD  +  5  =  0 


B 


08) 


so  that  the  desired  result  may  be  achieved  by  setting 

P  =  -2ek2  (19) 

Equations  (16)  and  (19)  suffice  for  the  determination  of  a ?  and 
ag  .  Using  these  and  the  definitions  contained  in  eqn.  (15)  and  in  [1], 
we  have 


a7  =  sg2bl 


a8  = 


ill 


2a2kB 


(20) 


where  the  parameters  b-j  ,  a-j  and  are  determined  as  in  [1]  and, 
as  previously  explained,  we  have  set  ac  =  1  . 

O 

The  approximate  dispersion  spectrum  is  now  completely  defined  and, 
from  (14),  the  dispersion  relation  is  found  to  be 
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(21) 


w2  ~  [-(0k^  +  yk2  +  <S  +  <$k2)  +  /D"]/2a 

J  J  A 


where 


0  =  (Sk^  +  yk2  +  6  +  $k2)2  -  4a(?k2  +  okx  +  0ky  +  pkx  ky  +  tIcx 

+  k®  +  ek2  kj)  (22) 

The  form  of  D  may  be  considerably  simplified,  however.  Making  use  of 
the  definitions  set  forth  above  and  in  [1],  we  may  show  after  some 
algebra  that 

<£2  -  4ot  5  0  ;  ~  2cte  =  0  ;  y<!>  -  2ap  =  0  ;  6<p  -  2aa  *0  (23) 

Thus  the  expression  for  the  discriminant  D  becomes 

D  »  (3kJ  +  yk2  +  6)2  -  4«Uk2  +  ek*  +  kj)  (24) 

and  we  see  that  the  discriminant  is  not  a  function  of  the  wave  number 
kx;  indeed,  it  reduces  to  a  form  identical  to  that  obtained  in  the  one¬ 
dimensional  theory  [1].  We  now  note  from  (16)  and  the  last  of  (23)  that 


<j>  ,  q 
Za  '  6 


so  that  (21)  may  be  rewritten  as 


(25) 


2 

U) 


s2k2 
9  ^ 


+  f(V 


where  f(k  )  is  the  double-valued  function 


(26) 
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(27) 


fUy  =  [-(3kJ  +  yk*  +  5)  ±  ✓D  ] /2cl 

and  D  is  given  by  (24). 

Necessary  and  sufficient  conditions  were  determined  in  [1]  such 
that  f(ky)  _>  0  for  all  real  values  of  k^.  This  requirement  was 
necessary  to  insure  real  frequencies  and,  indeed,  supplied  the  ration¬ 
ale  for  a  part  of  the  matching  procedure.  It  can  be  seen  from  (26) 
that,  for  f(ky)  >_  0,  the  frequency  u>  will  always  be  real -valued 
in  the  present  case.  For  constant  values  of  k^,  the  spectral  lines 
given  by  (26)  describe  hyperbolae  in  the  u  -  k^  plane.  In  the 
special  case  ky  =  0,  the  lowest  branch  degenerates  to  the  straight 

line  u  =  s„k  . 

g  x 

NUMERICAL  RESULTS  AMO  DISCUSSION 

As  a  numerical  example,  the  parameters  corresponding  to  the  exact 
solution  presented  in  [4]  were  used  to  evaluate  the  approximate  disper¬ 
sion  equation  (26).  Here  we  have  introduced  the  nondimensional  quanti¬ 
ties  a  =  2hcu/ircT  ,  c  =  2hK  ,/ir,  and  n  =  2hx  /*,  where  h  is  the 
i  x  y 

hal f- thickness  of  the  layers  with  unprimed  constants  and  Cy  *  /p/p  . 
Figure  2  shows  both  exact  and  approximate  lines  of  constant  ;  on  the 
dispersion  surface  over  the  first  two  Brillouin  zones.  In  the  first 
Brillouin  zone,  the  agreement  is  quite  good  through  ?  =  0.20.  In  the 
second  Brillouin  zone,  the  agreement  is  quite  good  as  well  for  g  =  0,  but 
worsens  as  c  increases.  These  tendencies  are  also  evidenced  in 
Figure  3,  which  shows  the  lines  of  constant  n  bounding  the  first  two 
Brillouin  zones. 
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Figures  4a,  b  and  c  show  the  real  and  imaginary  part  of  n  on  the 
imaginary  branch  connecting  the  first  two  Brillouin  zones.  Note  that 
the  exact  value  of  Re(n)  is  constant  on  those  branches.  The  agreement 
for  c  =  0  is  quite  good,  but  becomes  worse  with  increasing  c  ,  though 
still  adequate  through  5  =  .10  . 

On  the  basis  of  these  results,  it  appears  that  the  approximate 
theory  outlined  in  this  paper  yields  excellent  results  for  values  of  c 
and  n  up  to  about  one  Brillouin  zone  width,  with  somewhat  less  satis¬ 
factory  agreement  over  the  second  Brillouin  zone.  This  degree  of 
agreement,  however,  is  equal  or  better  than  the  results  of  other  exist¬ 
ing  approximate  theories.  In  addition,  the  present  theory  contains  the 
stopping  band  between  the  first  two  Brillouin  zones,  a  feature  present 
in  few  other  approximate  theories.  For  these  reasons,  the  two-dimensional 
effective  dispersion  theory  discussed  herein  is  thought  to  represent  a 
significant  advance  in  the  construction  of  approximate  theories  for  wave 
propagation  in  periodically  layered  bodies. 

It  should  be  noted,  however,  that  the  two-dimensional  effective 
dispersion  theory  is  still  subject  to  the  shortcomings  and  limitations 
of  the  one-dimensional  theory.  Principally,  these  are  a  somewhat  re¬ 
stricted  validity  due  to  the  requirement  that  the  frequencies  resulting 
from  the  approximate  dispersion  equation  be  real-valued,  and  the  lack 
of  a  uniqueness  proof.  Both  these  points  are  discussed  at  some  length 
in  [1]. 

We  close  by  briefly  discussing  the  prospects  for  improving  and 
extending  the  present  theory.  It  seems  likely  that  considerably  better 
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agreement  between  exact  and  approximate  solutions  can  be  obtained  by 
including  additional  terms  in  the  strain  energy  density  function  pro¬ 
portional  to  higher-order  derivatives  of  wQ  and  t|»  ,  and  by  intro¬ 
ducing  a  more  elaborate  matching  scheme.  This  improvement,  however, 
will  come  at  the  expense  of  greatly  increased  algebraic  complication. 
Additionally,  it  will  become  increasingly  difficult  to  satisfy  the 
requirement  that  the  approximate  dispersion  equation  yield  only  real¬ 
valued  frequencies.  Another  possible  improvement  is  the  extension  of 
the  theory  to  the  plane  strain  case.  Here  the  difficulties  are  due  to 
the  nature  of  the  exact  spectrum  for  an  unbounded  body,  which  has  been 
shown,  [5],  [6],  to  be  not  only  double-valued,  but  extremely  complicated 
as  well.  It  is  far  from  clear  as  to  how  a  successful  modeling  of  this 
very  complex  spectrum  ought  to  proceed. 
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FIGURE  CAPTIONS 

1.  Geometry  of  layered  solid 

2.  Curves  of  constant  ? 

3.  Curves  of  constant  n  at  Brillouin  zone  ends 

4.  Imaginary  branch  between  first  two  Brillouin  zones 

a)  c  =  0 

b)  c  =  .10 

c)  ?  =  .20 
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